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Several cosmic ray experiments have recently measured excesses in electrons and positrons, relative 
to standard backgrounds, for energies from ~ 10 GeV - 1 TeV. These excesses could be merely due to 
new astrophysical sources, but a more exciting explanation in which the electrons and positrons are 
dark matter annihilation or decay products is also consistent. Fortunately, the upcoming Fermi-LAT 
diffuse gamma ray measurements can further test these models, since the electrons and positrons 
produce gamma rays in their interactions in the interstellar medium. Although the dark matter 
gamma ray signal consistent with the local electron and positron measurements should be quite 
large, as we review, there are substantial uncertainties in the modeling of diffuse backgrounds and, 
additionally, experimental uncertainties that make it difficult to claim a dark matter discovery. In 
this paper, we introduce an alternative method for understanding the diffuse gamma ray spectrum 
in which we take the intensity ratio in each energy bin of two different regions of the sky, thereby 
canceling common systematic uncertainties. For many spectra, this ratio fits well to a power law 
with a single break in energy. The two measured exponent indices are a robust discriminant between 
candidate models, and we demonstrate that dark matter annihilation scenarios can predict index 
values that require "extreme" parameters for background-only explanations. 

PACS numbers: 



I. INTRODUCTION 

Since the introduction of dark matter into the scientific 
consciousness by Zwicky over 85 years ago, its nature has 
remained elusive. Though the existence of dark matter 
(DM) has been verified by numerous observations, we are 
still ignorant about its most basic properties, such as its 
mass and its interactions with standard model particles. 
A number of direct and indirect detection experiments 
are currently underway in an effort to change this, and 
several cosmic ray experiments have released interesting 
results over the past year. If these results are interpreted 
as signals of dark matter, then they lead to some sur- 
prising conclusions about this basic component of our 
universe. 

PAMELA [J, 0, i] has measured a pronounced in- 
crease in the positron fraction above 10 GeV. The sec- 
ondary positron spectrum, which originates from cos- 
mic ray (CR) interactions with the interstellar gas, is 
softer than the primary electron spectrum, so in the ab- 
sence of additional primary sources, the positron frac- 
tion was expected to drop with energy. The upturn 
in the positron fraction suggests the existence of a new 
primary source of positrons with energies in the range 
~ 10 — 100 GeV. Astrophysical sources such as pulsars 
i, & ! 0, H 1 EI El El offer one explanation for the 
PAMELA excess, but a dark matter explanation is an 
exciting possibility. 

For the total e + + e~ spectrum, ATIC [l3| found a 
broad excess in the range ~ 300 — 700 GeV with a sharp 
cutoff at the high energy end of the excess; PPB-BETS 
[l4| measured a similar excess. More recently, the Fermi 
collaboration [l5[ released its measurement of the e + +e~ 
flux over a very similar energy range. The Fermi spec- 
trum does not have the peak and subsequent sharp cutoff 



observed by ATIC, though it does exhibit a hardening 
around ~ 80 GeV. Between 100 and 500 GeV it is well 
described by a power law with an index that is harder 
than what is predicted by conventional theories of elec- 
tron propagation in the Galaxy. Additionally, the Fermi 
data suggests a break at the highest energies. HESS 
[l6l . ll~7T | has confirmed a break in the spectrum around 
~ 1 TeV. The results of all of these experiments are con- 
sistent with a new, primary source contributing to the 
local fluxes of electrons and positrons in the energy range 
10 - 1000 GeV. As with the PAMELA results, these may 
be a sign of dark matter. 

Regardless of the details of the production method, a 
population of high energy electrons and positrons will 
give rise to gamma rays with energies up to several hun- 
dred GeV through the process of inverse Compton scat- 
tering (ICS) off the low energy photons in the interstel- 
lar radiation field (ISRF). The Fermi Gamma-ray Space 
Telescope experiment is currently making all-sky mea- 
surements of gamma rays from tens of MeV to more than 
300 GeV with precise energy resolution and unprece- 
dented angular resolution. Fermi is therefore expected 
to detect the gamma ray products of the surprising new 
population of e^. The challenge will be to distinguish 
the gamma ray signal of the newly discovered high en- 
ergy e from the conventional background signal, so that 
the nature of the source can be ascertained. 

Several theories have been proposed to provide consis- 
tent explanations of both the PAMELA and Fermi results 
in terms of dark matter annihilation or decay [l8l . [l9l . I20I 

el m m, m, m m m m m m m, m m . while the 

details of the various models differ, explanations of the 
high energy signals as arising from DM annihilation 
have one thing in common: they require the annihila- 



2 



tion rate to be 2-3 orders of magnitude larger than that 
of a conventional thermal relic, i.e., the effective ther- 
mally averaged annihilation cross section today needs to 
be ~100 - 1000 times larger than 3 x 10" 26 cm 3 s _1 , the 
expected thermally averaged annihilation cross section 
for a WIMP with mass m x ~ 500 GeV. While this poses 
theoretical challenges, it is good news on the experimen- 
tal side, as it makes observing a dark matter signal at 
other indirect detection experiments more likely. Specifi- 
cally, if the additional high energy electrons and positrons 
observed at PAMELA and Fermi are a DM annihilation 
signal, then the observed large annihilation rate makes 
the prospect of measuring the ICS signal promising. 
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Many anal yse s of the DM gamma ray signal have been 
done [51. 135LI3R [37l. l38l. l39l. liol l-4ll| focusing on searches, 
for example, in dark matter subhalos, dwarf galaxies, and 
in the extragalactic diffuse emission. The galactic diffuse 
emission has long been the focus of dark matter gamma 
ray searches [H, |H, 0, Hi, Hf|. In spite of the large 
dark matter signal rate associated with the PAMELA 
and Fermi results, given an observed diffuse gamma ray 
spectrum, it is not straightforward to show that it re- 
quires a dark matter contribution. Unfortunately, theo- 
retical uncertainties in the background gamma ray spec- 
tra, along with the statistical and systematic uncertain- 
ties of the experiment, ensure that this is a challenging 
endeavor. Given the situation, it is useful to consider 
other approaches to analyzing the data. In this paper, we 
investigate an alternative method that removes some of 
the common systematic uncertainties that plague inter- 
pretations of gamma ray spectra. The idea is to consider 
the spectra for two different sky regions and take their 
intensity ratio in each energy bin, whereby systematic un- 
certainties that are energy dependent will cancel. More- 
over, fitting this ratio's energy dependence to a power law 
with a single break gives two measured exponent indices 
that can be robustly compared with different candidate 
models. As we demonstrate, for gamma ray spectra due 
to some dark matter models which fit the Fermi electron 
spectrum and early gamma ray data, the power law fit 
gives values for the indices that are hard to reproduce 
with background alone, unless rather extreme conditions 
are satisfied by the background. Since these indices are 
well measured and free from systematics, a signal consis- 
tent with the dark matter model can be more rigorously 
argued as evidence for its observation. 



The rest of the paper is organized as follows: in Sec- 
tion [TT1 we give a review of the different processes and 
components of the diffuse gamma ray spectrum. In Sec- 
tion mi we state our calculation methods, discuss our 
signal and backgrounds, and then describe in detail the 
ratio method. Readers who wish to skip these prelimi- 
naries should start at Section IIII Dl In Section IIVI we 
present our main results and in Section El we conclude. 
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FIG. 1: The components of the gamma ray sky: n 
(solid), ICS (dashed), bremsstrahlung (dotted), extragalactic 
isotropic (dot-dashed). The point source and DM contribu- 
tions are not shown. 



II. GAMMA RAY SPECTRUM 

The gamma ray sky as measured by Fermi is made up 
of several components: the 7r°, inverse Compton scat- 
tering and bremsstrahlung components coming from in- 
teractions of cosmic rays with the interstellar medium 
(ISM), as well as the point source and extragalactic 
isotropic components, as shown in Figure [TJ Addi- 
tionally, dark matter is expected to contribute, cither 
through direct production of gamma rays as it annihilates 
or decays, or through the interactions of its annihilation 
(decay) products with the ISM. 

The spectra of the various components can vary with 
the direction of observation, as the spectral shape and 
magnitude of a given component may be position depen- 
dent. The contribution from the prompt decay of neutral 
pions (ir° — > 77) produced in the p-p collisions of cosmic 
ray (CR) protons with the interstellar gas traces the dis- 
tribution of CR protons and the gas distribution, which 
is peaked toward z = 0. 

CR electrons, which include both the primaries ac- 
celerated by supernova shocks and the secondaries pro- 
duced in the p-p collisions of CR's with the interstellar 
gas, propagate through the Galaxy losing energy through 
through synchrotron radiation in their interactions with 
the Galactic magnetic field as well as inverse Compton 
scattering off of the photons of the Galactic interstel- 
lar radiation field (ISRF). Consequently, the distribution 
of ICS gamma rays is related to the distribution of the 
CR electron and the ISRF densities. The interstellar ra- 
diation field is composed of three sources: the low en- 
ergy, isotropic CMB photons, the optical photons due to 
starlight, which are concentrated near the Galactic Plane 
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(GP), and the infrared (IR) photons emitted by dust that 
has absorbed optical photons. The IR and optical com- 
ponents of ICS are thus very strongly peaked toward the 
GP where the starlight and dust are located. 

There is a diffuse isotropic component of gamma rays 
from unresolved sources such as active galactic nuclei and 
from diffuse emission processes outside of the Galaxy. A 
preliminary measurement of this cxtragalactic emission 
has been made by the Fermi collaboration and is com- 
patible with a power law index of —2.45 between 200 
McV and 50 GeV [13] as displayed in Figure □ Addi- 
tionally, there is the contribution from standard gamma 
ray point sources, such as pulsars and blazars, which are 
more concentrated near the Galactic Plane. 

Models of dark matter with high energy electrons 
and positrons as annihilation or decay products pre- 
dict a diffuse contribution to the gamma ray signal from 
bremsstrahlung and ICS [HI, H^, [Hfl, [5l|, [52|, [53|, [54| and a 
prompt contribution from final state radiation (internal 
bremsstrahlung) [H, M, M, E3, El, M, M, EO, 11 H • 
The prompt component has the same distribution as the 
square of the dark matter density in models of DM anni- 
hilation (or as the dark matter density itself in models of 
DM decay). Electrons and positrons lose energy through 
ICS with an energy loss rate that is proportional to the 
square of their energy. Therefore, the characteristic dis- 
tance they diffuse before losing most of their energy is 
much smaller for the highest energy electrons than for 
lower energy electrons. As a result, the highest energy 
ICS photons have a distribution that is similar to that 
of the prompt DM gammas, while the lower energy ICS 
photons exhibit a more diffuse distribution. 



cussion of the ISRF as modeled in GALPROP see [68|. 
Additionally, the atomic (HI), molecular (H2), and ion- 
ized (HII) hydrogen gas distributions are needed to com- 
pute the gamma rays from bremsstrahlung and 7r°-decay. 
The realization of these distributions in GALPROP is ex- 
plained in detail in [69]. 

For our calculations, we assume cylindrical symmetry, 
as well as mirror symmetry with respect to the Galactic 
Plane. Our standard diffusion zone has a radius of 20 kpc 
and a height of ±4 kpc. We take as our standard diffusion 
constant K{E) = 5.30 x 10 28 (£/4 GeV) 43 cm 2 s" 1 . 
These choices of diffusion parameters give values of the 
ratios B/C, sub-Fe/Fe, and 10 Bc/ 9 Be that are in agree- 
ment with the local data. 

We use a Galactic magnetic field parametrized as fol- 
lows: 



(B 2 ) 1/2 = B exp 



(r-Rp) jfj 
tb z b 



(1) 



where tb — 10 kpc is the scale radius, zb — 2 kpc is the 
scale height, and B = 5.0 pG is the local value of the 
field. 



B. DM Parameters 

We take the cusped Einasto profile as defined in Mer- 
ritt et al [7(| as our "benchmark" dark matter density 
profile. Additionally, we consider a cored Isothermal pro- 
file as defined in 64[ as an example of a cored DM density 
profile. These are defined as follows: 



III. ANALYSIS 
A. Galprop 

We use the GALPROP code (version v50p) of 
Moskalcnko and Strong to numerically propagate cosmic 
rays and calculate the resulting gamma ray production in 
the Galaxy [82j. In the spatial dynamics it includes dif- 
fusion resulting from cosmic rays scattering on magneto- 
hydrodynamic waves, while in momentum space, energy 
losses from ionization, bremsstrahlung, inverse Comp- 
ton scattering and synchrotron radiation are included 

Using the most recent measurements of the source 
abundances for the primary CR species, including nu- 
clei, electrons and gamma rays, GALPROP propagates 
the primary CR's through the Galaxy, iteratively com- 
putes the resulting spallation source functions for all 
species, and then propagates the full (primary + sec- 
ondary) source function for each species until a converg- 
ing result is obtained, assuming free escape of the parti- 
cles as the spatial boundary condition [H, H3] • 

The calculation of the ICS gamma rays requires knowl- 
edge of the interstellar radiation field. For a full dis- 



p(r) = po exp 



P(r) = Po- 



T->< '■ 



Einasto 



(2) 



Cored Isothermal 



(3) 



Rq = 8.5 kpc is the solar distance from the Galactic cen- 
ter, .13 < a < .22 is the parameter that defines the cuspi- 
ness of the Einasto profile (we use 0.17), r_2 = 25 kpc is 
the radius at which the logarithmic slope of the Einasto 
profile is —2, and r c = 2.8 kpc is the core radius for the 
cored Isothermal profile. We take po = 0.3 GeV cm~ 3 as 
the local value of the DM mass density. 

We define the boost factor (BF) of the dark matter 
signal as the ratio of the thermally averaged annihila- 
tion cross section needed to fit a set of data, (<7 Qnrl |w|) fu, 
to 3 x 10~ 26 cm 3 s _1 , the expected thermally averaged 
annihilation cross section for a WIMP with mass m x ~ 
500 GeV, 



BF = 



{Vann\v\)fit 

3 x 10- 26 cm 3 s 



(4) 



4 



C. Difficulties with Seeing the DM signal at Fermi 

As mentioned earlier, the dark matter gamma ray sig- 
nal detectable at Fermi consists of several components; 
the diffuse ICS and brcmsstrahlung signals from the in- 
teractions of the electrons and positrons produced in an- 
nihilations with the interstellar radiation field and the in- 
terstellar gas, respectively, and the prompt signals from 
internal bremsstrahhmg off of charged annihilation prod- 
ucts, i.e final state radiation (FSR), and from the decay 
of neutral pions produced in the annihilations [H| . 

For the dark matter scenarios we consider here, the 
diffuse brcmsstrahlung signal is 1-2 (or more) orders of 
magnitude smaller than the diffuse ICS signal, so we ne- 
glect it in our calculation. The FSR signal becomes rel- 
evant only at energies of order the DM mass (above the 
currently stated Fermi range of 10 MeV-300 GeV), and 
then only for certain annihilation modes, DM masses, 
and boost factors. Therefore, we neglect the contribu- 
tion of final state radiation to the dark matter gamma ray 
spectrum. Furthermore, we don't consider the effects of 
a prompt 7r°-decay component on the DM signal. Such 
a signal has a pronounced bump at high energies that 
extends well into the Fermi energy range and would be 
easily detectable for the relatively large branching ratios 
that would generically occur. However, the PAMELA 
collaboration's measurement of the antiproton fraction 
[7H places limits on the branching ratio to hadrons which 
severely constrain the flux of 7r°'s coming from hadronic 
annihilation channels. Fortunately, such suppressed or 
nonexistent ir° production can occur naturally in many 
dark matter models. Additionally, the total gamma ray 
flux as already measured by EGRET [zl, Izl and HESS 
[zH l7 5| p laces constraints on aprompt component of any 

type gi, m, m mi, m, m, izg. in h g ht f a n of this, 

for the dark matter signal, we henceforth consider the 
gamma rays produced from inverse Compton scattering 
of produced in DM annihilations. 

The Fermi e++e~ data suggest that M DM > 1000 GeV 
and that the spectra of from DM annihilation are 
fairly soft. For example, direct annihilation into elec- 
trons gives too hard a spectrum, but channels with softer 
electrons, like \X ~~ > an d XX ~* <t>4> followed by 

4> — > l + l~ , give good fits. Given the "large" annihila- 
tion cross sections these channels need to fit the Fermi 
and PAMELA data, BF > 200, one might expect the 
DM gamma ray signal to be easily detected at Fermi. 
However, the gamma ray signal may not be as large as 
the "local" signals suggest. The boost factor is highly 
position dependent, as it depends both on the extent of 
substructure in a particular region of the Galaxy and on 
the velocity dispersion of the dark matter particles, and 
the line-of-sight average may be significantly lower than 
the local value. 

Even with a large DM signal, it may prove difficult 
to disentangle the DM signal from the flux due to more 
mundane sources. The Galactic Center region, where the 



DM gamma ray flux is peaked due to the p 2 DM (r) depen- 
dence of the signal, has large contamination by gamma 
rays from unidentified point sources. A measurement of 
the DM component in this region depends on an accurate 
determination of the point source contribution. 

Perhaps the largest barrier to the detection of a 
DM signal is the uncertainty in the diffuse gamma ray 
background, the 7r°-decay gammas and the ICS and 
bremsstrahlung gammas from primary CR e~ and sec- 
ondary CR e ± . The diffuse background has uncertain- 
ties that arise from our incomplete understanding of cos- 
mic ray propagation, and from our inability to accu- 
rately measure the spectrum of cosmic rays throughout 
the Galaxy, the distribution of gasses in the interstellar 
medium, the Galactic magnetic field, and the distribu- 
tion of radiation energy density, which is closely tied to 
the distribution of dust in the Galaxy. Within these un- 
certainties, in many regions of the gamma ray sky, espe- 
cially those less prone to source contamination, it is often 
possible to take a spectrum that looks to contain a sig- 
nificant dark matter component and fit it to background 
alone. Thus, it is very difficult to use a single interesting 
gamma ray spectrum as evidence for dark matter. 

D. The Details 

We take the background gamma ray flux to be the 
sum of the contributions from 7r° decay, ICS, and 
bremsstrahlung and henceforth call this the ISM back- 
ground. The ICS and bremsstrahlung are themselves 
composed of gammas arising from the interactions of the 
primary CR e~ with the ISM, we call these the "leptonic" 
component of the ISM background, and of gammas aris- 
ing from the interactions of the secondary CR e with 
the ISM. Since these secondary e ± are produced in the 
interactions of cosmic ray protons with the ISM, we call 
the sum of these with the 7r° component the "hadronic" 
component of the ISM background. The leptonic com- 
ponent of the ISM background depends on the spectrum 
of primary electrons, while the hadronic component de- 
pends on the spectrum of primary protons. 

While there arc many DM annihilation modes that give 
good fits to the Fermi and PAMELA data, we use the 
annihilation mode through a light mediator <fi into muons, 
XX "~ *• 4>4> ~ > 2/! + 2/i~, as a template dark matter signal. 
Modes producing much harder e^- spectra don't fit the 
Fermi e + + e~ data as well, and softer modes produce 
similar spectra in the ~ 100 — 800 GeV range (though 
for different boost factors), as they must in order to fit 
the Fermi data. Since the populations of DM are 
similar, the dark matter ICS signals are therefore similar 
in the Fermi energy range. 

We calculate the gamma ray spectrum for several sce- 
narios, some with and some without a dark matter com- 
ponent. The ISM hadronic and leptonic components are 
generated by primary proton and electron spectra that 
are consistent with local data. We will describe in more 
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detail these local fits below. We fit the total flux, ISM 
+ DM ICS, to the gamma ray data recently presented 
by the Fermi collaboration [47|. The data are given 
for three regions, (0° < £ < 360°, 10° < \b\ < 20°), 
(0° < I < 360°, 20° < \b\ < 60°), and (0° < t < 360°, 
60° < which henceforth will be referred to as the 
Intermediate Latitudes 1 (IL1) region, the Intermedi- 
ate Latitudes 2 (IL2) region, and the Polar region, re- 
spectively, or collectively as the "regions". We assume 
that the extragalactic isotropic component and the point 
source component are well measured by the Fermi col- 
laboration, so we fit only to the Galactic diffuse flux as 
defined in the presentation of the Fermi "regions" data. 

We constrain the primary proton spectra to be in 
agreement with the local proton measurements by AMS- 
01 [z3|, BESS and IMAX [z| above 1 GeV. Ac- 

cordingly, we take the best-fit power law index below 9 
GeV to be —1.40. Above 9 GeV, we allow the power law 
index to vary such that the resulting spectrum does not 
exceed the local data by more than 3a. The indices for 
the fits to data of the hardest and softest allowed spectra 
are —2.00 and —3.50, respectively, while a very good fit is 
achieved for ~ —2.65. With the propagation and energy 
loss parameters we use in our calculations, the proton 
spectrum softens upon propagation by ~ —0.35 above 10 
GeV. Therefore, we take the range of acceptable power 
law indices for the injection spectrum of primary protons 
to be -1.65 to -3.15 above 9 GeV. 

Similarly, we constrain the spectrum of primary elec- 
trons in the range ~ 10 — 1000 GeV using the recent 
Fermi data Fitting to the data between ~ 73 GeV 
and ~ 380 GeV, we find that the range of indices —2.83 
to —3.20 gives spectra that are consistent with the data 
such that no point exceeds any Fermi data point by more 
than 3cr [Hj]. We stress that these constraints arise from 
fitting only the primary electrons to the e + + e~ data. 
Inclusion of a dark matter component of e^- in the flux 
allows for considerable softening of the allowed primary 
electron spectrum. However, we find that our conclusions 
are unchanged by considering a softer primary electron 
spectrum, so we take —3.20 to be the softest allowed 
power law index. Additionally, we note that we have not 
included the fluxes of secondary electrons and positrons 
in our fits, which are 1-2 orders of magnitude smaller 
than the flux of primary e~'s between 10 GeV and 100 
GeV. Including these would result in a small change of 
the spectral range. The electron spectrum softens upon 
propagation by ~ —0.65 in the energy range of interest 
for our choice of propagation and energy loss parameters, 
so we take the range of acceptable power law indices for 
the injection spectrum of electrons to be —2.00 to —2.55. 
All power law indices quoted for the electron and proton 
spectra will be those for the unpropagated spectra. 

With these choices, we believe we are being conserva- 
tive and are taking into account all of the realistic region 
of the parameter space of power law indices for both pri- 
mary protons and primary electrons, and perhaps some 
of the unrealistic region. While we don't think that the 
local spectrum of primary electrons is well-described by 
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FIG. 2: Plots for the contributions to the box region (10° < 
\i\ < 20°, 10° < \b\ < 20°) from dark matter (solid), lep- 
tons (dashed), and hadrons (dotted) with arbitrary normal- 
izations. Above 10 GeV, from top to bottom the index for 
the primary electrons and protons is —2.0, —2.5 and —3.0. 

the power law E~ 2 - s above 100 GeV, we do think that 
it is possible for the electron spectrum to be described 
as such in some region of the Galaxy, particularly near a 
CR source of electrons. Thus, we consider a wide range 
of hadronic and leptonic ISM spectra that are theoret- 
ically and experimentally motivated. See, Figure [2] for 
example spectra. 

As mentioned earlier, there are numerous sources of 
uncertainty in the diffuse background. For example, a 
change in the overall normalization of the magnitude of 
the Galactic magnetic field would increase the power to 
synchrotron energy losses, thus reducing the ICS signal 
everywhere, while a change in the parametrization of the 
magnetic field would ultimately lead to changes in the 
spatial distribution of the ICS. Similarly, changes in the 
optical and IR energy densities would result in a different 
ICS spectrum. As a final example, an increase in the HI 
density in a particular region would lead to a correspond- 
ing increase in the tt° component of the gamma rays. We 
don't attempt to address all of these issues here. Rather 
we expect that the large range of spectra we consider for 
the primary electrons and protons serves as a proxy for 
many of the uncertainties one expects from other sources. 



E. Taking the Ratio 

As our discussion on the gamma ray sources has made 
clear, it is challenging to conclusively determine that a 
measured gamma ray spectrum requires a dark matter 
contribution. Not only are there uncertainties in the 
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FIG. 3: The ISM background (dotted) and total, ISM back- 
ground + DM ICS, (solid) gamma ray signals for the Inter- 
mediate Latitudes 1 region (0° < I < 360°, 10° < \b\ < 20°) 
(upper) and Intermediate Latitudes 2 region (0° < £ < 360°, 
20° < \b\ < 60°) (lower). The plots are shown for comparison 
of the spectral shapes. 



background parameters, but experimental issues, in par- 
ticular energy normalization and systematic uncertain- 
ties, further cloud the issue. In this regard, it is useful to 
consider alternative analysis methods that can overcome 
some of these difficulties. 

The method we explore is to take the ratio in each 
energy bin of the gamma ray intensity for two different 
regions of the sky. This approach has several advantages. 
One, as displayed in Figure [TJ the standard backgrounds 
typically are dominated by a single contribution with a 
characteristic shape. For example, ir° gammas are ex- 
pected to dominate in the 1-10 GeV range in many sky 
regions. Thus, for energies where the two regions have 
the same process generating a majority of the gamma 
ray flux, the shapes look similar (sec Figure [3]) and give 
a roughly constant ratio; this in turn accentuates the en- 
ergy ranges where the gamma ray contributions are tran- 
sitioning from one source to another, making new sources 
like dark matter more prominent. Two, the systematic 
uncertainties that are common to a given bin, like the 
energy dependence of the acceptance, will cancel in the 
ratio. On the other hand, the overall energy normaliza- 
tion will still shift this ratio plot to the left or right, but 
since we choose to fit the ratio to a power law with a 
single break, this only shifts the break and not the ex- 
ponents of the power law. Similarly, the normalization 
of the power law is sensitive to overall factors that scale 
the different regions like any asymmetric exposure or ac- 
ceptance. Still, we expect that for a specific choice of the 
regions for the ratio, a power law fit with a break gives 
two measurable exponents that are a robust quantitative 



handle on the source of the gamma rays. 

As a specific implementation of this method, we take 
the average intensity of the "boxes" (10° < \£\ < 20°, 
10° < \b\ < 20°) and of the "strips" (\£\ > 20°, 10° < 
\b\ < 20°) as our two regions. These regions are well 
separated from the galactic center, making point source 
contamination more manageable, and for our analysis we 
assume that this has been properly subtracted. If point 
sources in regions closer to the Galactic Center, for e.g. 
5° < |6| < 10°, can also be reliably subtracted, choos- 
ing these closer regions would have an even larger dark 
matter signal, making the effect more pronounced. 

IV. RESULTS 

Our procedure is to take into account the statistical 
Poisson errors from one year of running at Fermi. We do 
not take into account systematic errors, as we argue they 
should largely cancel in the ratio. We bin in energies 
from 100 McV to 1.887 TeV using logarithmic binning, 
where the center of each bin is 1.2 times the previous. 
The acceptance is taken to be as given in [HJ multiplied 
by a relative change in selection (see slide 7 in [47[ )■ Note 
that at high energies, the expected observed gamma rays 
can be less than 1, leading to errors larger than 100%, 
which occurs above ~ 300-700 GeV for the box and strip 
regions. However, such bins are not crucial to the fits, so 
we leave them in to illustrate the behavior of the func- 
tion at high energy. A more realistic approach would be 
to look for an optimized binning of these high energy re- 
gions. Finally, we do not take into account any effects 
due to energy resolution. 

For most of the models, we first fit to the IL1, IL2, and 
Polar gamma ray data, marginalizing over a common set 
of indices and normalizations for the proton and electron 
primaries in these regions. However, for the dark matter 
models, we fix the electron primary index to be —2.45 
since this is used to determine the boost factor necessary 
for fitting the Fermi electron spectrum, but still allow its 
normalization to vary. Next, with the fitted parameters 
we predict the ratio of the intensity in the boxes over that 
in the strips and fit to a power law with a single break. 
Since the boxes have more dark matter signal than the 
strips, the dark matter contribution causes the ratio to 
increase at high energies. We include in the ratios the 
diffuse isotropic component of the gamma rays as given 
in |4?| . although it is usually unimportant. 

Our reference dark matter model, Model Rdm, is a 3 
TeV dark matter particle x annihilating in the channel 
XX ~~ * 00 2/i + 2/i _ with a boost factor of 990. This 
model gives a very good fit to the Fermi electron spec- 
tra, and the spectra in the IL1, IL2 and Polar regions 
can be fit well by varying the proton's normalization and 
index in addition to the primary electron normalization. 
In particular, this model's spectra in the IL1 region is 
compared to the best background fit to the regions in 
Figure |U In Figure [51 we plot the predicted ratio for this 
model with and without (Model Ro) the dark matter con- 
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FIG. 4: The preliminary Fermi-LAT gamma ray spectra in 
the Intermediate Latitudes 1 region (hatched band) as given 
in [47l ]. Three model fits are plotted: 1) the reference dark 
matter model Rdm (solid), 2) the best background only fit to 
the three regions BG reg ions (dashed), and 3) a hard electron 
spectrum to mimic the dark matter ratio behavior BG ra tio 
(dot-dashed). 



tribution. As can be seen from the plot, both ratios are 
compatible with a power law with a single break as shown 
by the solid lines. For both models, the ratio is roughly 
constant up to 1 GeV. On the other hand, from 1-200 
GeV the dark matter ratio is clearly increasing, while for 
the model without dark matter the ratio stays relatively 
constant. In Table HI we list the results of the fits to the 
gamma ray spectra in the IL1, IL2, and Polar regions 
and the power law fit to the ratio. For the indices of the 
power law, we've included the la error bars by marginal- 
izing the fit over the other parameters of the power law. 
The high energy index for the reference DM model is 
0.195^ 005, which is well separated from the value for the 
model with no dark matter, 0.049^ g 14 . 

As another comparison, we consider the best fit to the 
regions with background alone, Model BG rog i ns, which 
has primary proton and electron indices of —2.38 and 
—2.33, respectively. The ratio is shown in the bottom 
plot of Figure O The values for the low energy and high 
energy indices turn out to be similar to Model R , the 
reference model with dark matter turned off, as listed in 
Table Q] Incidentally, it is also worth mentioning that a 
background model that is specially designed to fit very 
well the Fermi electron spectra gives a high energy index 
of 0.038. Thus, it seems generic that the best-fit back- 
grounds to either the Fermi electrons or the Fermi IL1, 
IL2, and Polar gamma rays do not generate a high energy 
index that is as large as that for dark matter. 

We investigate some variations in the dark matter sce- 
nario to determine the effects on the values of the power 



FIG. 5: Ratio plot for Model Rdm (top), 3 TeV dark matter 
with BF — 990, and Model Rn (bottom), same background 
with no dark matter contribution, as a function of energy in 
GeV. The fit to a power law with a single break is shown in 
solid. 



law fit. First, we reduce the boost of the dark matter 
to (1/2,1/3,1/4) which changes the high energy index 
values to (0.130,0.107,0.094) with only a small change 
to the low energy index. Second, we change the mass 
of the dark matter particle to 1 TeV and 5 TeV, tak- 
ing, respectively, boosts of 200 and 2000 to give a rea- 
sonable fit to the Fermi electron spectrum. In these 
cases, we find high energy index values of 0.113 and 
0.240, with low indices of -0.035 and -0.014. Third, 
we consider a cored Isothermal dark matter profile with 
BF = 530, and we find a high energy index of 0.154. 
This is not surprising; the cored Isothermal profile is 
less cuspy than our standard Einasto profile, and so pro- 
duces a smaller gamma ray flux near the Galactic Cen- 
ter. Finally, we alter our diffusion parameters from our 
standard values. We consider a diffusion constant of 
K{E) = 1.29 x 10 28 (£/4 GeV) 40 cm 2 s" 1 with a dif- 
fusion zone of \z\ < 1, and also a diffusion constant of 
K(E) = 7.00 x 10 28 (£/4 GeV) 43 cm 2 s" 1 with a diffu- 
sion zone of \z\ < 6 [85] . Here we find the high energy 
index to be 0.035 and 0.199, respectively. This makes 
physical sense, because for the latitudes we are looking 
at, \z\ > 1.5 kpc at r = 0, hence the smaller diffusion 
zone is suppressing the electrons from dark matter anni- 
hilations that would contribute to the gamma ray flux. 
As all of these variations demonstrate, there is some un- 
certainty in the high energy index value. However, to be 
optimistic we continue to use the reference scenario as 
the target signal we are searching for. 

It is important to determine if it is possible to vary 
the background to consistently fit the gamma ray spec- 
tra in the regions while getting similar index values for 
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Model 


Description 


Proton 
Index 


Electron 
Index 




ffibins 


Low Energy 
Index 


High Energy 
Index 


-^break 

(GeV) 


Rdm 


Reference model 
with dark matter 


-2.47 


-2.45 


0.23 


0.33 


-0.0171;°°$ 


0.195l;° ° 


1.15 


Ro 


Reference model 
without dark matter 


-2.47 


-2.45 


0.65 


0.09 


-0.046!;°°,? 


0.0491;°^ 


0.93 


BGregions 


Best BG fit to IL1, IL2, 
and Polar regions 


-2.38 


-2.33 


0.11 


0.07 


-0.034t;gE 


0.049l;° J° 


0.93 


-° ^regions 


Best BG fit to regions with 
modified norms in boxes 


-2.38 


-2.33 


0.11 


0.98 


-0.1331;°°* 


0.137t;°» 


0.83 


BGratio 


Hard electron spectrum 
fit to ratio 


-3.00 


-2.00 


0.54 


0.67 


-0-0001;°°^ 


0.204l;°°J 


0.91 


BG rat ; 


Slightly softer electron 
spectrum fit to ratio 


-2.52 


-2.20 


0.13 


0.91 


-o.o89±:SSS 


0.179l;°°l 


0.83 



TABLE I: The models with their indices for the primary proton and electron spectra (unpropagated), along with their fit 
parameters to the IL1, IL2, and Polar regions, and for the ratio plot, the parameters for a fit to a power law with a single 
break. 
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FIG. 6: Ratio plot for Model BGregions (bottom), the best 
background fit to the IL1, IL2, and Polar regions, and Model 
BGr Cgions (top), same as Model BGregions with a change in 
relative normalization of 1/2 and 2 for the primary proton 
and electron spectra in the box region. The fit to a power law 
with a single break is shown in solid. 



the power law fit. This is one of the benefits of the ra- 
tio method, as the power law indices give a quantitative 
set of numbers to compare to that are not limited by 
systematics. For our simplified survey into the possible 
background variations, we not only vary the primary pro- 
ton and electron spectra, but also allow for a change in 
their normalization in the boxes relative to the normal- 
izations used in the IL1, IL2, Polar, and strip regions. 
Such normalization changes could be mocked up by vari- 



ations in Galactic quantities in the boxes relative to the 
strips. In fact, to first order, changing the gas density and 
the strength of the IR component of the ISRF would, re- 
spectively, mimic a change in the hadronic and leptonic 
normalizations. Based on matching the reference model's 
high energy index alone, we find that these relative nor- 
malizations must be allowed to be as large as a factor of 
2. 

Our study of these variations shows that there is a ten- 
sion in the background achieving a good fit to both the 
regions and the power law fit, while also obtaining both a 
large high energy index and small low energy index. The 
reason for this is due to the particular energy dependen- 
cies of the different contributions as shown in Figure [2] 
Notice that the dark matter ICS contribution (solid) is 
growing as a function of energy with a broad plateau 
from a few GeV to 100 GeV, whereas the backgrounds 
have a much different behavior. The hadronic contribu- 
tion (dotted) has a "camel hump" plateau from about 200 
MeV to 1 GeV with a high energy tail that falls faster for 
softer values of the primary proton index. For the leptons 
(dashed) , the contribution is a smoothly falling function. 
In general, the region fits favor the hadronic component 
to dominate at its plateau, with some leptonic contri- 
bution at higher energies. Since the leptonic component 
does not grow as a function of energy like the dark matter 
component, to get a large high energy index for the power 
law fit to the ratio, one needs the lepton component in 
the boxes to be enhanced relative to the hadronic com- 
ponent. Due to the shape of the leptonic contribution, 
it will then also give a sizeable contribution to the low 
energy region below the peak of the hadronic component. 
Thus, the box region's spectra in increasing energy will 
go from being dominated by leptons, to being a mixture 
of leptons and hadrons around the hadronic plateau, to 
again being dominated by leptons. Since there are three 
regions of interest, the fit to a power law with a single 



9 



10 




1 - 



0.1 1 10 100 1000 

Energy (GeV) 



FIG. 7: Ratio plots for Model BG ra tio (top), a background 
with indices of —3.00 and —2.00 for primary protons and elec- 
trons, respectively, with an increase in relative normalization 
of 1.6 for the primary electrons in the box region, and for 
Model BGratio (bottom), a background with indices of —2.52 
and —2.20 for primary protons and electrons, with a change 
in relative normalization of 0.6 and 1.7 in the box region. The 
fit to a power law with a single break is shown in solid. 



break will have a poorer % 2 than the case of dark matter 
which is dominated by hadrons at low energy and dark 
matter at high energy. In the cases where the leptons 
strongly dominate at low energies, the value of the ra- 
tio will actually fall as a function of energy, leading to a 
noticeable negative low energy index. 

To illustrate these conclusions, we look at some specific 
cases. We consider a case, Model BG rcgions , in which we 
take the best fit background to the regions and set the 
hadronic normalization in the boxes to be 1/2 and the 
leptonic normalization in the boxes to be 2 times that 
of the strip region. The ratio is shown in the top plot of 
Figure [6] It is noticeably falling at low energies and rising 
at high energies. As seen in Table [H both indices are well 
separated from reference model Rdm , and the \ 2 f° r the 
power law fit is markedly worse. However, the low and 
high energy index can be brought into agreement with the 
reference model by a judicious choice of indices for the 
protons and electrons. Taking softer protons and harder 
electrons brings out the high energy lepton contribution 
without as large of a correction at low energies. 

In Model BGratio we take the proton and electron pri- 
mary indices to be —3.00 and —2.00, respectively, corre- 
sponding to a very hard electron spectrum. This choice 
allows for a good balance between decent fits to the 
gamma ray data in the IL1, IL2, and Polar regions and 
a ratio that mimics the one for the DM reference model. 
For a lepton normalization in the box 1.6 times that in 
the strip, we get values for the low and high indices sim- 



ilar to the DM reference model, as seen in the top plot 
of Figure [7] A distinguishing feature of this model is 
that in the region plots, there is a dip in the spectra 
around 3 GeV, visible in FigureSl near the transition be- 
tween hadronic and leptonic contributions. In addition, 
the high energy part of the ratio plot has a distinct long 
plateau starting at 30 GeV, which is considerably differ- 
ent from the reference dark matter model which has a 
narrow plateau starting at about 200 GeV. As reflected 
by the % 2 fits, with further statistics it will be possible 
to test if either of these features are truly there. 

Finally, we note that the challenge of mimicking the 
dark matter signal increases if we consider a stricter con- 
straint on the electron index from the fit to the Fermi 
electron data. As our Model BGr atio , we take the elec- 
tron index to be —2.20 and the proton index to be —2.52 
to better fit the region plots. With a modification of the 
hadron and lepton normalizations of 0.6 and 1.7 in the 
box, as shown in the bottom line of Figure [7J we get a 
reasonably large value of the high energy index, but also 
a large negative low energy index. 

Although this exploration of the background variations 
has been quite simplistic, it illustrates that the power 
law fit to the ratio gives useful quantities for models to 
compare to. In our study, we have shown that it is nec- 
essary to take somewhat extreme parameters for the pri- 
mary proton and electron spectra to get a ratio that is 
similar to that for a scenario with a dark matter con- 
tribution. For one, the electron (proton) index typically 
needs to be on the hard (soft) end of the range of its 
allowed values. In addition, their normalizations in the 
box region have to be as much as a factor of 2 different 
from that inferred in the strip region. Although more 
systematic background studies should be performed, we 
believe that our general conclusions would still be valid. 
This is because for the energies we are interested in, the 
hadronic and leptonic gamma ray contributions are dom- 
inated by a single process. Reasonable variations of these 
background components will then largely be captured by 
our normalization and index changes. Therefore, we con- 
clude that there is, quite generically, tension in achieving 
a background model that fits well the gamma ray data 
for the IL1, IL2, and Polar regions and that has a large 
high energy index and small low energy index. Thus, if 
Fermi were to observe the ratio with a similar behavior 
as the dark matter reference model, this could be confi- 
dently argued as evidence for a dark matter signal in the 
gamma rays. 

V. CONCLUSIONS 

The future Fermi-LAT measurements of the diffuse 
gamma ray spectrum have the potential to resolve the 
degeneracies between standard and dark matter explana- 
tions of the cosmic electron and positron excesses. How- 
ever, background and experimental uncertainties make it 
difficult to claim a dark matter discovery even if the spec- 
trum is consistent with such a model. In this regard, it is 
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useful to have many different approaches to this analysis. 
In this paper, we suggest a new technique in which the 
intensity ratio is taken for gamma ray spectra in two dif- 
ferent sky regions. Some systematic uncertainties cancel 
in the ratio, and by fitting the ratio to a power law with 
a single break, we find that the two measured power law 
indices give a robust quantitative handle of the data. 

Operationally, we fit our models to the preliminary 
Fermi gamma ray spectra in the Intermediate Latitudes 
1, Intermediate Latitudes 2, and Polar regions. We then 
predict the gamma ray spectra in the box and strip re- 
gions that comprise part of the Intermediate Latitudes 
1 region. Computing the ratio of these two spectra, we 
show that the dark matter scenario has a relatively flat 
ratio at low energies and an increasing ratio at high ener- 
gies. The flatness of the ratio at low energies reflects the 
similarity of the hadronic contributions to the spectra in 
the two regions, while the increase in the ratio at high 
energies reflects the increase in the dark matter contribu- 
tion at high energies in the box region, which is closer to 
the center of the Galaxy. Thus, for the power law fit, the 
signal has a small low energy index and a large high en- 
ergy index. By comparison, background explanations of 
the Fermi electron and gamma ray data predict a much 
smaller high energy index. The background ICS signal, 
having a spatial dependence that is quite different from 
that of the DM ICS signal, has a much smaller enhance- 
ment in the box region relative to the strip region. 

We also investigate the types of background parame- 
ters leading to spectra that could mimic the dark matter 
behavior, and find that it requires quite extreme values. 
One, the electron and proton primary spectra must be, 
respectively, near the hardest and softest spectra allowed 
by local data. Two, the electron normalization in the box 
region must be enhanced relative to the proton normal- 
ization by a factor of ~ 2 — 4. In general, it is difficult for 
background-only models that fit well the IL1, IL2, and 
Polar regions' gamma ray data to have a ratio similar to 
the ratio obtained in the dark matter scenario. To ob- 
tain a high energy index for a background-only model, 
one must invoke a very hard leptonic component, which 
is marginally consistent with the local data and possi- 
bly inconsistent with expectations for source production. 
Given this ability to quantitatively compare models and 



robustly demonstrate the background's difficulty in im- 
itating the dark matter behavior, a Fermi observation 
consistent with the dark matter signal can be more con- 
fidently argued as evidence for dark matter. 

Further studies of this method should be performed 
to see if these preliminary results are generic. For the 
background, many different parameter variations can be 
studied beyond our simplified look at adjusting proton 
and electron primary indices and normalizations. For 
instance, pulsar explanations can be studied to deter- 
mine their predictions. For signals, dark matter models 
beyond our "reference" model can be analyzed, in par- 
ticular additional annihilation channels as well as dark 
matter decay scenarios. It is possible that the ratio anal- 
ysis will help make models with much smaller gamma ray 
signals detectable. Additionally, the method can also be 
further optimized. In this paper, we focused on the in- 
termediate latitude box and strip regions to compute the 
ratio; other choices may be preferred, while using several 
regions and taking many ratios could provide even more 
useful information. For example, if point source contri- 
butions can be reliably subtracted for spectra at lower 
latitudes, the ratio of these regions would have an even 
more pronounced dark matter signal and require more ex- 
treme backgrounds to explain them. On the other hand, 
ratios of two background-dominated regions may be use- 
ful in understanding the backgrounds, helping to reduce 
the background uncertainties. With future studies and 
applications along these lines, this new analysis technique 
may prove useful in explaining cosmic ray excesses and 
even in providing further evidence for dark matter. 
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